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ABSTRACT 

The complex-time formalism is developed in the framework of the path- 
integral formalism, to be used for analysis of the quantum tunneling phe- 
nomena. We show that subleading complex-time saddle-points do not ac- 
count for the right WKB result. Instead, we develop a reduction formula, 
which enables us to construct Green functions from simple components of 
the potential, for which saddle-point method is applicable. This method 
leads us to the valid WKB result, which incorporates imaginary-time in- 
stantons and bounces, as well as the real-time boundary conditions. 
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1. Introduction 



The imaginary-time path- integral method has been successful for the treatment 
of the quantum tunneling phenomena in quantum field theories. The existence of 
the solution of the field equation, instantons and bounces, allows us to apply semi- 
classical approximation valid in small coupling. 

In recent years, much effort has been made to develop this method further. 
One of the driving force was the hope of observing baryon and lepton number 
violation process in collision experiments in the TeV range. This process is as- 
sociated with the tunneling between topologically different vacua in the standard 
electrowcak model through the baryon and lepton number anomaly, 't Hooft'^' 
first noticed this possibility, but at the same time found that the probability is 
highly suppressed. Later it was found that the potential barrier that separated 
the vacua has a "pass" (the Sphaleron), at the height of few TeV. This implies 
that even if the tunneling at the bottom of the potential well is highly suppressed 
as was estimated by 't Hooft, it may not be suppressed once initial energy is of 
order of TeV. This expectation lead to series of works on developing the original 
instanton methods to be applicable to higher energy tunneling. This development 
has important implications in quantum theory in general, including nuclear fission 
problems as well as a number of tunneling problems in matter physics. 

The original dilute-instanton-gas approximation is valid at the ground state. 
This can be confirmed by comparing the result in quantum mechanics with that 
of the WKB calculation. As the initial energy goes up, instanton result becomes 
invalid. In quantum mechanics, it fails in the estimate of the energy splitting of the 
first (and higher) excited state in double well. In quantum field theory, the cross 
sections violate the unitarity.'^' One of the present authors (H.A.) and Kikuchi'"' 
found that the interactions between the instantons become relevant: It gives the 
right energy splitting in quantum mechanics and satisfies the unitarity bound for 
the cross section. 

As the energy yet goes up, the distances between instanton and anti-instanton 
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decreases. Thus the semi-classical configurations are far from the solution of the 
equations of motion. Consequently, the valley methods were adopted. The new 
valley method'*^' was developed for this purpose in mind and yields configurations 
that allows exact conversion of the smallest eigenvalue to collective coordinate. 
Khoze and Ringwald used a series of approximate valley configurations to es- 
timate the total tunneling cross section. The new valley method was recently 
apphed to the bounce problem by one of the authors and Wada. They solved for 
the bounce-valley and identified small bubble configurations valid for the finite 
energy decay. 

All these analyses were carried out in the imaginary-time formalism. Thus basic 
as well as practical questions related to analytic continuation remain. How does the 
boundary condition matches that of the valley configurations? In estimating many- 
point Green functions, how does the external fields enter (affect) the imaginary 
time configuration? [This is crucial for the validity of the path-deformations in 
the valley calculation, as well as bounce calculations.] These are some of the 
relevant questions. Some other points are elaborated by Boyanovsky, Willey and 
Holman!'*' They also note that in some cases the imaginary-time method leads us 
to apparently contradictory results. 

In quantum mechanics, the complex-time method has been studied by various 
authors!"* It was argued that it allows semi-classical approximation for tunneling 
phenomena and overcomes various problems associated with the pure-imaginary- 
time method!^' In it, one considers the analytic continuation of the time- integral in 
Fourier transform of the Feynman kernel. The existence of the tunneling leads to 
the existence of the complex saddle-points in the time-plane. It was claimed that 
semi-classical approximation is done by deforming the time-integral so that it goes 
through all such saddle-points.'^^' 

This gives one hope that calculations in field theory may be improved by using 

[17—19] 

the complex-time method. Indeed, Son and Rubakov adopted a complex-time 
method for the estimate of the cross section. They used periodic instanton-anti- 
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instanton solution in imaginary time, and the initial Minkowski state overlaps 
with the instanton as a coherent state, much like the Aoyama and Goldberg'^"' 
calculation of the Sphaleron cross section. 

[21] 

Closer look, however, reveals many unanswered questions. It is known that 
the complex-time plane is plagued with singularities and infinite number of saddle- 
points, which form lattice structure. Among saddle-points, there are two kinds, 
those that are solutions of the equation of motions (we shall call these "physical 
saddle-points" in this paper), and those that are not ("unphysical saddle-points"). 
How the path is deformed to avoid the singularities and to go through only the 
physical saddle-points is unknown. One simply assumes that all and only the 
physical saddle-points contribute. Even so, the weight of each saddle-point is a 
riddle. They are determined so that the result agrees with that of the WKB 
approximation. It is not known how it comes about from the path. 

In this paper, wc give the solution of these problems, based on a reduction 
formula in the number of the turning points of the path. In the next section, we 
give a brief overview of the saddle-point calculations in complex-time. There we 
elaborate some of the points briefly mentioned above. In section 3, we analyze 
the orthonormahty of the WKB wavefunctions and construct the Green function 
from them. The connection formula for the WKB wavefunction is discussed in the 
section 4, where it leads to the reduction formula for Green function. Expanding 
this formula, we obtain a series, which can be interpreted as a sum over the physical 
saddle-points. The weights and the phases of saddle-points arc determined from 
this formula. The last section gives summary and discussions. 
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2. Saddle-point method 



We consider one-dimensional quantum mechanics with action; 

'1 



dt 



- V{x) 



(2.1) 



The potential V{x) is assumed to be smooth enough to allow WKB approximation 
in asymptotic regions I (x ^ — oo) and II (x ^ oo). 

The finite-time Green function (Feynman kernel) is defined by the following in 
the Heisenberg representation; 



G{xi,xf,T) = {xj\e-'"'^\xi), (2.2) 

in terms of the hamiltonian of the system, H. The advanced and retarded resolvents 
are, respectively; 





G\xi,Xf-E)=t J dTe'^^-'^^'^G{xi,xf,T) = {xf\^—^^—j^\xi), 

— oo 

(2.3) 

oo 

G'^ixuxf-E) = -zj dTe'^''+'^^TG{xi,xf,T) = {xf\ ^ _ ^ 



\x. 



where real-positive 5 is introduced to guarantee the convergence of the integrals. 

The poles of these Green functions come from the bound states, and the cut from 
the continuous spectra. 

Let us first look at the saddle-point approximation for the path integral, 

G{xi,Xf;T)^ J Vxe'^. (2.4) 

x{Q)=Xi 
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One obtains the following from the saddle-point method; 



G{xi,Xf;T) = J2 



Xcl 



2ti Xd{0)xa{T) 
i -d'^Sci/dT^ 



.-1/2 



(2.5) 



where x^i is the solution of the classical equation of motion and Sd is the action 
evaluated at the classical solution. This approximation is valid as long as both 
Xi and Xf are far from the turning points where classical velocities vanish. Sub- 
stituting (2.5) into the retarded resolvent in (2.3), one arrives at the following 
expression; 



oo 

G^{xi,xf,E) ^ -i J dTj2 



2n XcmxdjT) 
i -d^Sd/dT^ 



-1/2 



i{ET+Sci) 



(2.6) 



In the above, we absorbed the infinitesimal imaginary convergence factor (iS) in 
E. 

In calculating (2.6), one hopes to apply the saddle-point approximation to the 
T-integral. The saddle-point condition (the stable-phase condition, to be more 
exact) is then. 



E 



9^cl 1-2,1// ^ 

-g^ = ^Xd +V{Xd) 



(2.7) 



In other words, the solution of the equation of motion Xci{t) has to have the energy 
E. 

Let us consider a situation depicted in Fig.l. Both the initial point Xi and 
the final point xj lies in the forbidden region. There is no real-time solution of 
the equation of motion. Thus no semi-classical approximation is possible for the 
real-time formalism. 

If one assumes that the expression (2.5) gives correct analytic continuation 
of the Green function in the complex T-plane when complex-T solutions of the 
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equation of motion exist, one could go to the complex T-plane and apply the 
saddle-point method. One then arrives at the following expression; 

G(xi,Xf;E) = -iJ2J2\^cimcl(Ts)\-'/' w(Ts) p(Ts) e^^-(^^), (2.8) 

% Xcl 

where Tg is given by (2.7). Here Wci{Ts) is the WKB phase, w(Ts) is weight and 
p{Ts) is phase for the corresponding saddle-point. The original contour is deformed 
to pass the saddle-points by a series of steepest descent paths from the origin T — 
to the positive infinity on the real axes (Re T = oo, Im T = 0). The weight of a 
saddle-point, w{xci), is determined by how the contour crosses the saddle-point: If 
the contour does not cross the saddle-point, the corresponding weight is zero. If 
the contour crosses it along the steepest descent direction, the weight is 1. In case 
the contour reaches a sub-dominant saddle-point from the direction orthogonal to 
the steepest descent and leaves it along the steepest descent direction, the weight is 
1/2 (as in the calculation of the false vacuum decay). The phase of a saddle-point, 
p{Ts) comes from the square root in Eq.(2.6). 

Carlitz and Nicole applied this method to linear potential, quadratic well and 
quadratic barrier, all of which are exactly solvable. Using the exact expression 
for the Green function (2.2), they found that the saddle-points explained above, 
namely the ones associated with the complex-time solution, indeed lead to saddle- 
points in the T-plain [we shall call these "physical saddle-points"]. They, however 
found more: there were saddle-points that are not associated with any solutions of 
equations of motion (unphysical ones), and also singularities due to the periodici- 
ties. Deforming the T-integration path to avoid singularities, they found that the 
path indeed goes through only the physical saddle-points. As a result, the weights 
and phases in (2.8) are determined, which lead to the correct WKB result. 

In the case of a double-well potential, the deformation of the integration con- 
tour is not specified, which is understandable in view of the fact that the analytic 
structure of (2.2) in complex T-plain is fairy complicated. The claim made in the 
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literatures is simply that the contour passes all saddle-points that correspond to 
classical trajectories with weights that are extracted from the above calculation. 

The problem lies in the fact that we do not know the contour which passes 
all physical saddle-points that are distributed in complex-time plane. Because of 
this, we cannot determine the weights of the saddle-points. In the case where Xi 
and X f are in the forbidden region on the right side of the double- well potential as 
depicted in Fig.l, saddle-points in complex T-plane are given by 



Ts = ±T{b2, Xi) ± T{b2,Xf) + lT{ai, bi)+mT{bi, 03) + nT{a2, 62) 

1,171,11 — 0, ±1, ±2, ... , 

where 



(2.9) 



y I 

T{x,y) = 2 / (2.10) 



X 



Notice that T(x, y) is pure imaginary if {x, y) is in a forbidden region. These saddle- 
points Ts arc shown as solid circles in Fig. 2. The unphysical saddle-points arc also 
shown as open circles in the same figure. We know of no contour which passes 
only the physical ones and avoids the unphysical ones. Thus so far the complex- 
time method has not been worked out for the double-well and more complicated 
potentials. 

3. Green function and the WKB wavefunctions 

We shall examine the Green functions by using the complete orthonormal set 
of the eigenfunctions {|'0n)} of H. For example, the retarded Green function in 
(2.3) is written as 

G {xi, xf, T)^l^ -^T^53a;^' (3-1) 

n 

where the sum over the complete set is actually made of the integration over the 
continuum spectrum and sum over the discrete bound states. 
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In order to study of the Green functions in the asymptotic region, we use the 
WKB approximation for the wavefunctions ip{x) = {x\ip). For a state with the 
eigenvalue A in the continuum spectrum, the second-order WKB approximation 
yields the following; 

il^xix) = ^ X \ ' (3.2) 

^/^^ I {Cxe" + Dxe' ^(^')'^"') for x G II, 



where p{x) = -^2(A — V{x)). We choose the limits xi and a;2 of the integrations to 
be the nearest turning points for definiteness. Among the coefficients A, B, C, and 
D, only two are independent. Their inter-relation is linear due to the superposition 



[22] 

principle and can be written as follows; 




(3.3) 



The 2x2 matrix S{X) is determined by the shape of the potential V{x) in the 
intermediate region between xi and X2- The flux conservation law is satisfied for 
the Schrodinger equation; 

d . „ . i f dip* 

Substituting the WKB expression (3.2) into the above, we find that 

j = l.^\A\'_\B\^) = ^{\C\^-\D\'). (3.5) 
This results in the following relation between the matrix elements of -S'(A); 

|^ll|'-|^2l|'=l, |5l2|'-|-S22|' = l, 

(3.6) 

'-'ll'-'12 — '-'21'-'22 — U. 

As a result, the matrix S{X) is parametrized by three real functions a{X),(3{X), 
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and p(A) as follows; 



e*" cosh p e*^ sinh p 
S{\) =1 I . (3.7) 

e sinh p e cosh p 



For a given eigenvalue A, there are two independent eigenf unctions. In order to 
construct them, we look at the inner- product of the eigenfunction ipx{x) with coef- 
ficients {Ax, B\, C\, D\) and another eigenfunction ip\'{x) with {Axr, By, C\i, Dy). 
Integration in the asymptotic regions determines the coefficient of the delta func- 
tion 6{X — A') completely, as the delta function can come only from the infinite 
integrations. [In order to calculate the finite terms we need to solve the Schrodinger 
equation completely. But this is not necessary for the current purpose.] In fact, 
integration in the region II yields. 



oo oo 

If dx 



J dx 'iIjI{x)'iI;x'{x) 



X2 Xl 



27r j wWxWTx) 

Xl 

X I (7*(7^,e*-C(p(*')-p'(a;'))da;' ^ ^^j^^^-ill^{jp(x')-p\x'))dx' (3.8) 



As we are interested in the singularity of the above for A = A', we expand the 
exponent using the following for A ~ A', 



V{x')-AA^^^^- (3.9) 



By using a new coordinate y defined by dy — dx/p{x), we find the contribution 
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from the first two terms to the delta function; 

CO CXD 

a;2 2/2 

= \{ciCx + DlDx) 5{\-\') + ... 

(3.10) 

where we have neglected all the finite terms. Doing the similar calculation for the 
region I, we find 

oo 

j dxrx(x)Mx) = ^ (a^Aa + BlBx + ClCx + DlD^) 5{X - A'). (3.11) 

— oo 

Using the above result, we choose our orthonormal eigenfunctions for a given A as 
the following; 



(„ia \ 
0, 1, -e^(«-^)tanhpj, 

e^(«+^)tanhp, 1, 0, j . 



(3.12) 



[The above corresponds to two incoming states. There are of course other choices, 
such as stationary states, but all of those reads to the same result for the Feynman 
kernel.] 

Let us first look at the case when both xi and Xf are in an allowed region in 
II. Prom (3.12), we find that 

j:4\xf)4>ix,)^ \ 

itii 27r y/p{Xi)p{Xf) 

Jxi > _q\Jx2^^ Jx2^^ ;g«(.a tanh p + (c.c.) 

(3.13) 

The complex conjugate part can be understood as the same function below the cut 
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on the real- axis of the complex A-plane (see Fig. 3). This is guaranteed by the fact 
that as A moves below the cut, the role of the coefficients A and B is exchanged 
and likewise for C and D, and as a result, the phases a and (3 change their signs, 
while p does not. Therefore the sum over the continuous spectra in the Feynman 
kernel (3.1) is written as the following A-integration; 



The integration contour C in the complex A-plane is as in Fig.3. In (3.14), we have 
chosen the integrand so that it converges for |A| oo, which we need later. As for 
the discreet spectrum, it is determined by the absence of the diverging behaviour of 
the wavefunction, which translates into the condition that for C = 0, A = 0. From 
(3.7), we thus find that sinhp = when A is equal to a discreet eigenvalue of H. 
This allows us to write the sum over the discreet spectra as pole integrations of the 
second term in (3.14) depicted in Fig.3. Connecting all the contours and closing 
it at the infinity, we find that the whole contour enclose the pole a,t X — E + iS. 
Thus we end up with the expression 





(3.14) 



{xi,Xf;E) 



I 



\/p{xi)p{xf) 




(3.15) 



where we defined the "reflection coefficient" R by the following. 



R = i^ =ie*("-^) tanh p. 



(3.16) 
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We note that this expression (3.15) is vahd for a general value of E with proper 
analytic continuation of the coefficients. Therefore, if any end-points are in a 
forbidden region, simple analytic continuation of (3.15) is appropriate. We shall 
use this result in the next section. 

4. Reduction formula and its expansion 

In this section, we derive a formula for evaluating the Green function in a 
system with an arbitrary potential to which we can apply the WKB approximation. 
We consider a general potential which has arbitrary number of wells, and show that 
the Green function is given by summing up all contributions of classical paths. 

Let us consider the case where Xi and Xf lie in the forbidden region on the 
same sides of the wells. The potential which has n wells is depicted in Fig. 4. By 
the analytical continuation of (3.15) in E, we find the retarded resolvent in the 
following form; 

G^{xi,Xf;E) = -\p{xi)p{xf)\-^/h-^' (e^f + tRne-^A (4.1) 

where 

Aij= J dx\p{x)\. (4.2) 

bn 

In order to specify the fact that this expression is for n wells, we attach the subscript 
(n) to the coefficients hereafter. 

Due to the existence of the intermediate region (bn-i,an), the matrix 5**^"^ 
which connects the regions I and II, can be written in terms of the matrices S^'^~^^ 
that connects regions I and {bn-i, an) and S for {bn-i, an) and II. If we apply the 
linear WKB connection formula for the latter region, which we can obtain from 
the saddle-point method for a linear potential, we obtain the following; 

5(")(B) = S(»-)(B)(^"^^"" 7) 
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X 



e 



l[Wn+%) ] 1 1 



2 



( ^^ I \e-^^-^cosWn e-^"-i sin \ 



where Wn and A„_i are defined by, 

6n 



iyn = y" (4.4) 



an 



A„_i= y da;|p(a;)|. (4.5) 

b„-i 

Therefore, we find the relations; 

S'i^ = -^g-'^e-^-i cos Wn - 4r'^e^"-i sin Wn (4.6) 

= -Sg-'^e-^-i sin + 25g-')e^"-i cos Wn . (4.7) 
Prom (3.16), (4.6) and (4.7), we find that Rn can be written as the foUowing; 

1 - sin W„ - ^e~^^"-ii?n_i cosl-Vn 

Rn^ ^^TT — 4.8 

2 cosl^n - 5e-2An-ii?„_isinWn ^ ^ 

We find it most convenient to rewrite this to the foUowing two expressions; 

11-/? p2iW^n 



i?^ = f . (4.10) 

1 + ii?n-ie-2^"-i ^ ^ 



1 -|^_ie-2A„-i 

Defined in this manner, the function Rn-i corresponds to the reflection amphtude 
at the turning point bn-i- 
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In order to see the correspondence between these expressions and the saddle- 
point method, let us expand (4.9) in an infinite series as follows; 



iRn^\+ [-iRn) e^^^" + {-iRn?{-i) e^^^" + ... . (4.11) 

The convergence of this series is guaranteed by the imphcit factor i5. The retarded 
resolvent is then, 

G^(a;„x-;;i?) = -|p(a;,)p(a;;)|-V2 

^ ,-(A,-A,)^|^^j^_^^^j e2^^"+(-z^n)'(-0e^''^" + -}e-(^^+^^) ■ 

(4.12) 

This is the expression to be compared to the saddle-point expression (2.6). If we 
choose T to be negative imaginary, — ir, the factor in the exponent is, 

i{ET + Sci) ^ Er - Se ^ - J dr (^^^ = - J dx\p{x)\ (4.13) 

Therefore, the two factors A^j is equal to the above quantity for the path from 
Xij to the turning point bn- The first term in (4.12) corresponds the contribution 
of the pure- imaginary-time path that starts from Xi and reaches x f directly. The 
factor 6"^^*+^^ ) in the rest of the terms is for the path from Xi to the turning point 
bfi and then from bnto xj. The expansion of iRn is understood as contributions of 
the paths that oscillate in the allowed region {un, bn)- Various factors have unique 
interpretations as factors coming from the turning points. This is most conveniently 
depicted in Fig. 5. Similarly to (4.11), the expression (4.10) is expanded as the 
following, 

-iR^ = -i + (ii?„_i) e-2^"-i + {iRn-ifil) e-^^"-' + - ■ (4.14) 
The corresponding diagrams are illustrated in Fig.6. 
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We have so far derived the expressions that results from the WKB approxima- 
tion for the n-th well. This procedure can be applied recursively for the rest of the 
wells. At the end, we are left with Ri, which is given by 

iRi = }-tanWi='- + i-i) e-2*^i + {-if e'^^' + ... (4.15) 

where W\ is defined by (4.4) by replacing n with 1. Therefore we find that Ri is 
expressed as a sum of the contribution of the classical paths which evolve in the 
first well. 

Combining (4.12), (4.14), and (4.15), we find all complex-time paths are in- 
cluded in the expression of the resolvent G^. Therefore, the resulting Green func- 
tion in a n-well potential is given by 

G'^ixi^xf^E) = -\p{xi)p{xf)\-^'^Y. /(^d)e*^^"""^'''^- (4.16) 

where f{xci) is determined by the number of refiections and transmissions which 
the classical path contains, and W{xi,Xf; E) is determined by what wells the path 
crosses and what barriers it tunnels through. The rules for calculating f{xci) is 
the following: When a path has a reflection in the allowed region, it obtains —i. 
In the case of a reflection in the forbidden region, it gets |. 

Let us next examine the case when Xi, Xf are on opposite sides of the wells. 
Applying the analysis similar to the one in the section 3, we find that the Green 
function is expressed as follows; 

G''{xi,Xf;E) = -\p{xi)p{xf)\-'/^Tne-^^-^f (4.17) 

where Tn is the (analytically-continued) transmission amplitude 

'-'22 

Just as the previous case, we can express Tn in terms of T^-i as in the following; 

T„ = I 

4i ^^e-^r.-ismWn + 2S^'^ ^^e^"-i cosM^n 
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rp -An-l+iWn 

J-n—l^ 



1 + ii?n-ie-2^"-i + e2iW^- (1 - ii?^_ie-2An-i) 

In the above, we have used (3.16), (4.7) and (4.10). Let us show that T„ consists 
of all contributions of classical paths. The second factor in (4.19) is expanded in 
the following way, 

^e^^- + {-i){-iRn)e'''^- + {-if{-iRnf e^'^- + ... , (4.20) 



1 + Rne^i^^ 

which corresponds to the diagrams in Fig. 7 (a). The first factor in (4.19) is ex- 
panded as. 



= e-^"- + (J)(ii?,_i) e-3^"- + (J)2(ii?,_i)2 e-5^"- + ... , 



(4.21) 

which is illustrated in Fig. 7 (b). As the previous case, repeating this procedure, 
we reach Ti, which is given by 



2cosW^i 

^ e'^' + {-if e^'^' + + ... . (4.22) 

Therefore we conclude that consists of all contributions of classical paths. 

The above analysis can be applied to other cases with various locations of Xi 
and Xf. Thus we verify that thus constructed resolvent can be always interpreted 
as sum over the complex-time classical paths. 

This formalism is valid for a metastable system and scattering process (Fig.8) 
as well as for a stationary system. It is apparent from this derivation that this 
complex-time method reproduces the result of the WKB approximation with the 
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linear connection formula. The potential we consider is given in Fig. 8(a). The 
turning points are labeled in the figure by ai, bi, and 02- First we consider the 
case when Xi is in the well and Xf is outside of the barrier. The analysis similar to 
the above yields the Green function as in the following; 

G\x„.r,E) = -i\p(xM^,)\-'l' \^^^^,,^^^^^^,^_,^^^ . (4.23) 
where R is given by 

1 - ie-2A^ 



—iR — —i 



A' 



1 + le-2Ai 



(4.24) 



= -i + le-2^^ + (l)3e-4^^ + .... 



From (4.23) and (4.24), we find again that the resolvent is equal to the sum over the 
physical saddle-points. We find that the poles of the Green function are determined 
by the following; 

1 + e^W^ + 1 (1 - e2^^^) = 0. (4.25) 

Let us solve this equation iteratively. Then we obtain 

W,{E)^{n+^^7:-\e-^^K (4.26) 

From (4.26), we find an imaginary part of the energy eigenvalue 

Im En = , (4.27) 

where T{En) is the period of the classical path between the turning points ai, 61 

fei 

/rlr 
, (4.28) 
V2 {E - V{x)) 

ai 

We find the decay rate of this metastable system. We note that the factor 1/2 in 
(4.27) comes from the weight 1/2" of the saddle-points with n reflections in the 
forbidden region. 
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Let us apply this method to a scattering process against a potential barrier in 
Fig.8. When xi and Xf are separated by the barrier, the Green function is given 

by 



G''{xi,Xf;E) = -ib(xi)p(x;)|-VW(^^+^/) 



k=l 



2(fc-l) 



(4.29) 



where we have used the WKB expression for the transmission coefficient T; 



-Ai 



1 + ie-2Ai 



(4.30) 



When Xi and x f are on the same side, the Green function is given by 



G''{xi,Xf;E) = -t\p{xi)p{xf)\-'/\''^^ 



-i\p{xMxf)\-'/\''^^ 



k=l 



(4.31) 



where we have used the expression for the reflection amplitude. 



—iR — —i 



1 - 4e 



1 + ie-2Ai 



(4.32) 



The coefficients T and R satisfy the unitarity condition 



lTr + \Rr = 1- 



(4.33) 



We again confirm the validity of the sum over the physical saddle-points. The 
weight 1/2" is crucial for (4.33). 
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5. Discussion 



In the paper, we have given the re-formulation of the complex-time method. 
Using the connection formula for the wavefunctions, we constructed the reduction 
formula in the number of the turning points for the Green function. This yields 
series expansions, which can be understood as sum over in the classical complex- 
time trajectory. This is understood as a sum over the physical saddle-points with 
specific weights and phases in the path-integral method. This shows the valid- 
ity of the method proposed before, in the context of the path-integral method. 
We confirmed that this method yields results identical with that of the WKB ap- 
proximation. Thus our construction gives solid basis for the starting point of the 
complex-time method. 

What happens if one tries to extend this method to field theories? Not much 
is known. Rubakov et.al. basically assumes that the behavior of the Green func- 
tions in asymptotic region (T — > oo) does not depend on the imaginary part of 
T and calculate with non-zero Im T so that the saddle-point approximation can 
be applied. However, even for Re T — > oo, the Green function does depend on 
the Im T, as we are allowed to have more configurations along the imaginary-time 
Im T increases. Thus rather careful consideration would be required. The 
complex-time method instead relies on the T integration as we have seen in this 
method. Thus somewhat different analysis may be required. One of them could be 
the combination of the valley method and the current complex-time method. The 
valley method can used to identify the imaginary-time tunneling paths, which are 
converted to collective coordinates. Thus, while the incoming particle are expressed 
in real-time expressions, the tunneling part may be obtained as the imaginary part 
of the complex-time development along the valley trajectory. There are rather 
interesting possibilities along this line, which should be further investigated. 



20 



REFERENCES 

1. G. 't Hooft, Phys. Rev. D14 (1976), 3432. 

2. N. S. Manton, Phys. Rev. D28 (1983), 2019; F. R. Klinkhamer and 
N. S. Manton, Phys. Rev. D30 (1984), 2212. 

3. A. Ringwald, Nud. Phys. B330 (1989), 1; O. Espinosa, Nud. Phys. B343 

(1990) , 310. 

4. H. Aoyama and H. Kikuchi, Phys. Lett. B247 (1990), 75, Phys. Rev. D43 

(1991) , 1999, Int. J. Mod. Phys. A7 (1992), 2741. 

5. A. V.Yung, Nud. Phys. B191 (1981), 47. 

6. H. Aoyama and H. Kikuchi, Nud. Phys. B369 (1992), 219. 

7. V.V. Khoze and A. Ringwald, Nud. Phys. B355 (1991), 351. 

8. H. Aoyama and S. Wada, Kyoto University preprint KUCP-0069 (1994). 

9. D. Boyanovsky, R. Willey and R. Holman, Nud. Phys. B376 (1992), 599. 

10. D. W. McLaughlin, J. Math. Phys. 13 (1972), 1099. 

11. I. Bender, D. Gomez, H. Rothe and K. Rothe, Nud. Phys. B136 (1978), 
259. 

12. S. Levit, J. W. Negele and Z. Paltiel, Phys. Rev. C22 (1980), 1979. 

13. A. Patrascioiu, Phys. Rev. D24 (1981), 496. 

14. A. Lapedes and E. Mottola, Nud. Phys. B203 (1982), 58. 

15. U. Weiss and W. Haeffner, Phys. Rev. D27 (1983), 2916. 

16. R. D. Carlitz and D. A. Nicole, Ann. Phys. 164 (1985), 411. 

17. S. Y. Khlebnikov, V. A. Rubakov and P. G. Tinyakov, Nud. Phys. B367 
(1991), 334. 

18. V. A. Rubakov, D. T. Son and P. G. Tinyakov, Phys. Lett. 287 (1992), 342. 

21 



19. D. T. Son and V. A. Rubakov, Nucl. Phys. B424 (1994), 55. 

20. H. Aoyama and H. Goldberg, Phys. Rev. Lett. 188B (1987), 506. 

21. H. Aoyama and T. Harano, Kyoto University preprint, KUCP-0072/KUNS- 
1302 (1994) 

22. H. Aoyama and M. Kobayashi, Prog. Theor. Phys. 64 (1980), 1045. 



22 



FIGURE CAPTIONS 



1) A potential V{X) with asymptotic regions I and II. 

2) Positions of saddle-points in complex T-plane. 

3) The complex-plane of the eigenvalue A. The cut on the real-axis is the con- 
tinuous spectra and the poles are the bound states. The integration contour 
C is for the continuous part of the Green function. The discreet part can be 
represented as integration around poles. 

4) A potential with n wells we consider. The turning points are denoted by 
and bi, for left and right of the i-th well, respectively {i — 1 ^ n). 

5) Diagrammatic representation of the expansion of iRn- The first term corre- 
sponds to the direct reflection at the turning point bn- Each path in the rest 
corresponds to the group of paths that have a number of oscillation in the 
n-th well. The information of the region left of is contained in —iRn- 

6) Diagrammatic representation of the expansion of —iRn- The each path cor- 
responds to a group of paths with different number of oscillation in the for- 
bidden region (6„_i,a^). 

7) Classical paths traversing in the n-th well and tunneling through the (n — 1)- 
th barrier. These are again expanded into groups of paths with different 
oscillations in either the allowed or forbidden regions. 

8) A metastable system and a scattering process against a potential barrier. 
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